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ABSTRACT 


The radiation mechanism of fast radio bursts (FRBs) has been extensively studied but still remains 
elusive. Coherent radiation is identified as a crucial component in the FRB mechanism, with charged 
bunches also playing a significant role under specific circumstances. In the present research, we propose 
a phenomenological model that draws upon the coherent curvature radiation framework and the mag- 
netized neutron star, taking into account the kinetic energy losses of outflow particles due to inverse 
Compton scattering (ICS) induced by soft photons within the magnetosphere. By integrating the ICS 
deceleration mechanism for particles, we hypothesize a potential compression effect on the particle 
number density within a magnetic tube/family, which could facilitate achieving the necessary size for 
coherent radiation in the radial direction. This mechanism might potentially enable the dynamic for- 
mation of bunches capable of emitting coherent curvature radiation along the curved magnetic field. 
Moreover, we examine the formation of bunches from an energy perspective. Our discussion suggests 
that within the given parameter space the formation of bunches is feasible. Finally, we apply this model 
to FRB 20190520B, one of the most active repeating FRBs discovered and monitored by FAST. Several 
observed phenomena are explained, including basic characteristics, frequency downward drifting, and 
bright spots within certain dynamic spectral ranges. 


Keywords: Radio transient sources (2008); Radio bursts (1339); Neutron stars(1108); Magnetars (992) 


1. INTRODUCTION 


'The systematic researches of distant fast radio bursts 
(FRBs), a class of bright millisecond radio transients, 
have started in 2007 (Lorimer et al. 2007). Throughout 
16 years of accelerated development, numerous signif- 
icant observational milestones have been accomplished 
(Spitler et al. 2016; Chatterjee et al. 2017; Bochenek 
et al. 2020; CHIME/FRB Collaboration et al. 2020), 
which are comprehensively reviewed in the literature 
(Cordes & Chatterjee 2019; Petroff et al. 2022). To date, 


Corresponding author: Chengmin Zhang (zhangcm@bao.ac.cn) 
and Di Li (dili@nao.cas.cn) 


over 700 FRBs (The CHIME/FRB Collaboration et al. 
2021) and 5,000 bursts! have been recorded, and several 
expansive datasets corresponding to individual sources 
have been published (Li et al. 2021; Niu et al. 2022; Xu 
et al. 2022). In addition, the abundant data also provide 
a foundation for statistical studies of FRBs (Feng et al. 
2022), such as FRB luminosity function (Luo et al. 2018; 
Hashimoto et al. 2022), FRB population analysis (Cui 
et al. 2021; Pleunis et al. 2021), FRB as cosmological 
probes (Macquart et al. 2020; James et al. 2022), etc. 
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Despite rigorous exploration, the precise origin and ra- 
diation mechanism of FRBs remain enigmatic. Regard- 
ing their origin, the magnetar model (Popov & Postnov 
2010) garners support from observations of the Galactic 
SGR 1935-2154 (Bochenek et al. 2020; CHIME/FRB 
Collaboration et al. 2020). However, no such bright 
bursts have been detected in subsequent observations 
(Lin et al. 2020), and comparable radio bursts have not 
been exhibited by other magnetars. As for their radia- 
tion mechanism, coherent radiation has emerged as the 
most feasible explanation for FRBs, given their exceed- 
ingly high brightness temperature (Lu & Kumar 2018; 
Lyutikov 2021), yet the exact physical processes are con- 
tinue to be a topic of considerable debate. 

Broadly speaking, the mechanism of FRBs can be 
divided into two categories: inside the magnetosphere 
(pulsar-like) and outside the magnetosphere (GRB-like), 
for detailed information, can see the model review 
(Lyubarsky 2021; Xiao et al. 2021; Zhang 2022). Among 
these, the model proposing coherent curvature radiation 
by bunches (Katz 2014; Cordes & Wasserman 2016; Ku- 
mar et al. 2017; Ghisellini & Locatelli 2018; Yang & 
Zhang 2018; Wang et al. 2020; Cooper & Wijers 2021) 
represents one of the foremost contenders for explain- 
ing FRBs. These models posit that FRBs stem from 
the curvature radiation generated by relativistic bunches 
traversing along the curved magnetic field lines. Ob- 
servational constraints stipulate that the size of these 
bunches (l) should be less than half of the wavelength 
(e.g., A/2 ~ 10cm for 1GHz) to ensure that the ra- 
diation signals from various charged particles within a 
single bunch reach the observer in approximately the 
same phase (Wang et al. 2020). This phenomenon is 
thus termed coherent curvature radiation. 

Nonetheless, these models confront two primary ob- 
stacles. The initial challenge pertains to the forma- 
tion of relativistic bunches (Melrose 2017), and the sec- 
ondary challenge relates to the influence of plasma ef- 
fects (Lyubarsky 2021; Beloborodov 2021). In address- 
ing the latter issue, Qu et al. (2022) demonstrated that 
bright FRBs could escape the magnetosphere across a 
vast parameter space and propagate outward. But, 
the first issue still lacks a satisfactory resolution. Al- 
though the hypothesis of a two-stream instability has 
been posited (Usov 1987; Melikidze et al. 2000), sim- 
ulation studies to corroborate this conjecture are ab- 
sent under extreme relativistic conditions (Samuelsson 
et al. 2010; Tokluoglu et al. 2018). More specifically, the 
formation problem is attributed to the possibility that 
Coulomb forces, electrostatic potential, or other inter- 
actions could cause dispersion among charged particles 


within the bunch, thereby inhibiting their congregation 
(Zhang 2022). 

In this paper, we aim to analyze the formation of dy- 
namic bunches through the inverse Compton scattering 
(ICS) mechanism, which precedes curvature radiation 
and leads to additional energy loss that cause compres- 
sion of the outflow stream. Building on this model, 
we intend to further investigate the observational fea- 
tures of FRB 20190520B. This FRB was detected by the 
Five-hundred-meter Aperture Spherical radio Telescope 
(FAST) (Li et al. 2018; Niu et al. 2022) and is distinct 
not only due to its association with a persistent radio 
source (PRS) (the other is FRB 20121102A), but also 
owing to its sustained activity since its discovery (Niu et 
al. 2023 in prep.). This allows for ongoing monitoring 
and the discovery of new features, and thus we have se- 
lected FRB 20190520B as a case study for our modeling 
and analysis. In Section 2, we establish the model from 
the perspective of the bunch compression process. In 
Section 3, we present simulations based on our model, 
applying this model to FRB 20190520B, and providing 
explanations for observed phenomena. Subsequently, in 
Section 4, we discuss other related implications and con- 
cerns. Finally, we summarize our conclusions in Section 
5. 


2. COMPRESSED BUNCH MODEL 


In this section, we will construct our model from a 
physical perspective. Initially, a general overview of the 
model will be provided, followed by the introduction of 
three model components. Subsequently, we present the 
derivation of the model and parameters, underpinned by 
plausible assumptions related to neutron stars (NSs). 


2.1. General illustration 


In our model, we consider the whole process of 
FRB emission, drawing upon the framework of coher- 
ent curvature radiation by bunches and the fundamen- 
tal physics of NS (Lorimer & Kramer 2012; Lyne & 
Graham-Smith 2012). Generally speaking, this kind of 
model should depict at least three components, as shown 
in Figure 1, which encompass the original energy source 
(trigger region), the formation of bunches (outflow re- 
gion), and the radiation ensuing from the movement 
of bunches or interaction with the surrounding medium 
(FRB region). 

Firstly, the trigger mechanisms for FRBs are generally 
associated with several instances of intense energy re- 
lease, such as crust cracking of the NS surface (Thomp- 
son & Duncan 2001; Beloborodov & Thompson 2007), 
sudden magnetic reconnection (Lyubarsky 2020), or the 
discharge of magnetic energy following glitches in the 
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Figure 1. Schematic diagram of one burst based on the compressed bunches model and coherent curvature radiation from a 
magnetized neutron star (NS). In sub-figure (a), the brown cylinders represent various bunches with different heights (hi & h2), 
Lorentz factors (Yo1 & Yo2), and curvature radii (p1 & p2). The red wavy lines signify the coherent curvature radiation emitting 
from the bunches. The yellow curves depict the magnetic tubes, and the blue dashed lines represent the line of sight (LoS). 
Moreover, the regions of trigger, outflow, and FRB emissions are also displayed. In sub-figure (b), the striped area circled by 
the orange dashed line represents the outflow stream that undergoes compression via inverse Compton scattering (ICS), leading 
to the formation of the dynamic bunch. The dark and light orange part represent the upstream and downstream of the outflow 
particles, respectively. T'he regions of pair production and ICS are also marked. 


NS (Wang et al. 2021). Thus, all these mechanisms can 
be considered as being magnetically powered. ‘These 
events, which we regard as feasible in the context of 
astrophysics, could result in an immense outflow of rel- 
ativistic charged particles. Theoretically, such events 
need to occur outside of the outflow magnetic field lines, 
such as beyond the open magnetic field region (Lu et al. 
2020). If this is not the case, the violent eruptions could 
disrupt the ordered magnetic field. Therefore, only a 
minor quantity of charges will enter the open field re- 
gion and stream out along the magnetic field following 
pair production, providing the seed charges for FRBs. 

Secondly, the formation of bunches poses a key ques- 
tion in this type of model (Kumar et al. 2017). In this 
context, we present a potential scenario for the forma- 
tion of bunches. In contrast to other bunched coher- 
ent curvature radiation mechanisms, we postulate that 
other radiation processes precede FRB emission (Lon- 
gair 2010), including pair production, bremsstrahlung, 
and inverse Compton scattering (ICS). These processes 
are situated closer to the NS surface (in the outflow re- 
gion) than the curvature radiation of FRB (FRB region), 
as depicted in Figure 1. Given the considerations of 
mean free path (acting distance) (Timokhin & Harding 
2019) and scattering cross-section, we argue that ICS is 
the predominant process (refer to Appendix B for fur- 
ther details). Therefore, before reaching the point of 
FRB emission, the outgoing relativistic charged parti- 
cles undergo energy loss due to ICS. 


Assuming an environment saturated with soft photons 
originating from the NS surface and pervading the in- 
ner magnetosphere, the energy loss rate of a relativistic 
particle becomes significant. Consequently, its kinetic 
energy experiences a considerable reduction even over 
a short period of ICS action. Simultaneously, for the 
upstream particles (the preceding portion of the out- 
flow stream, represented by the darker segment in sub- 
figure (b) of Figure 1), the duration of ICS is extended, 
resulting in a more pronounced energy loss compared 
to the downstream particles (the subsequent portion of 
the outflow stream, represented by the lighter segment). 
Macroscopically, the stream undergoes compression as 
it moves, tending towards unification for six dimen- 
sions (three spatial and three momentum dimensions). 
'This process may facilitates the dynamic formation of 
bunches. 

Finally, the compressed relativistic bunches traverse 
along the curved magnetic field, emitting coherent cur- 
vature radiation in the tangent direction with various 
heights. The frequencies of this radiation are correlated 
with the distinct radii of curvature and the Lorentz fac- 
tor (y) unique to each bunch. Among these, the num- 
ber density of the compressed bunch may be exceedingly 
high, thereby a magnetic tube/family should be consid- 
ered, which is related to the state of NS (Link 2003; Lan- 
der 2013). It is crucial to emphasize that the bunches are 
dynamically compressed during their movement and ra- 
diate at certain instants, without the need for long-term 
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stability. The previously mentioned content provides a 
phenomenological description of the compressed bunch 
model, and we will now proceed to outline the physical 
derivation of this model. 


2.2. Physical derivation 


Let us start with the first component: from the energy 
budget to the estimation of total particles. Assuming we 
interpret the central engine as a magnetized compact ob- 
ject, such as a magnetar (Duncan & Thompson 1992). 
Then, we can compute the approximate total magnetic 
energy within a specific volume as Erot ~ B?V,,/(87), 
where B, ~ 10!°G denotes the surface magnetic field, 
and Vm ~ 3 x 101? cm? signifies a volume that can com- 
pletely release its magnetic energy (refer to Appendix 
A for parameter estimations). The Frot should consti- 
tute the upper limit of a single FRB burst. The re- 
leased energy surges out of the trigger region as photons 
or electrons with an extraordinarily high Lorentz factor 
ye ~ 10° (Ruderman & Sutherland 1975; Medin & Lai 
2010). Then the total number of particles is, 


Niot = nBiVa/(81ymec)) 10? B3 15Vm 1370 (1) 


in which 7 is efficiency factor of energy conversion, Me 
is the mass of an electron, and c is the speed of light. 
If the conversion rate of this process is close to nuclear 
energy release, then the 7 ~ 107°. 

However, considering that the trigger region is out- 
side the outflow region, and not all outflow particles will 
be coherent in the FRB region, some constraint factors 
need to be further added. For example, ¢ is the parti- 
cle number rate from the trigger to the outflow region 
~ 10-?; also, € means the volume factor, which will 
be given later. Meanwhile, considering that the tube is 
not full of the whole outflow region, so another factor 
is & ~ 1078 (details see Appendix A). Therefore, the 
number of seed particles that can coherent is 


Nic = RCEB2VA/ (Bn^ymec?) 


2 
~ 1.5 x 107£ B2 15V. i815; (2) 


where € is given in Eq.8. When these ultra-relativistic 
particles flow out, their high Lorentz factor cannot be 
maintained because of the pair production process (Ru- 
derman & Sutherland 1975). The cascade is not the 
key question of this article, so we will only describe it 
phenomenologically. Generally, the change of the y is 
related to the cascade number, and the pair production 
becomes the main factor in the decrease of 7, such as 
in the condition of yg. However, as the y decreases to 
73 ~ 10? (Longair 2010), the effect of ICS becomes sig- 
nificant (see Appendix B). Thus, the cascade number is 


estimated as (Timokhin & Harding 2019) 


Neas = logo ^r [ ya) RI 10, (3) 


and the particle number after the pair production is 
Ne = 2Ncs Nic. Here, we have roughly completed the 
first part of the model that the estimation of particle 
number for the energy budget. 

Next, the second part focuses on the energy loss and 
compressed mechanism to dynamically form the rela- 
tivistic bunches. As we considered in Section 2.1, par- 
ticles have energy loss mechanisms during their out- 
flow process, and the energy variation can be expressed 
through changes in the y. If we believe that the up- 
stream particles lose energy for a longer time than the 
downstream particles, then the upstream part will ex- 
perience a greater energy loss and have a smaller y. In 
other words, there will be a speed difference between 
the upstream and downstream, and a distance change 
will occur after some time. Specifically, this manifests 
as a reduction in the distance difference between them, 
and the downstream will move towards the upstream, 
resulting in a compression effect, as illistrated in Figure 
1. Then the compressed length is written as 


8l ~ (v; — Uo) dt = (v; — Vo) (Yi — yo) Mec? / Eos, (4) 


where ôt is the compression time, v is the velocity of 
electrons, subscripts i and o represent the parameters 
at the beginning (in) and end (out) of the compression, 
respectively, and Pisas is the power of energy loss mech- 
anism. As analyzed above, the ICS is dominated during 
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Figure 2. Evolution path of the Lorentz factor (y) and 
compression length of the outflow particles with time (t) 
due to the inverse Compton scattering (ICS) compression 
effect. The blue, orange, and green solid lines represent the 
initial yrcs; with different values at the beginning of ICS 
energy loss (see left ordinate). The blue, orange, and green 
dashed lines represent the compression length (ôl) with dif- 
ferent values of yros: (see right ordinate). The red dash-dot 
line means y = 100. 


this process, which power for a single electron can be 


calculated as (Longair 2010) 


4 
Poic = 3°09 Upn; (5) 
where 8 = v/c ~ 1 under the relativistic condition, 


and Uy, = OTt /c is the internal energy of photons 
with a constant temperature distribution T ~ 10° K 
(Goldreich & Julian 1969) (surface temperature of a NS) 
and Stefan-Boltzmann constant os. If we consider that 
y varies with time, the Eq. 5 can be written as follows, 


Pec = Pioss = Mec’ Yo, (6) 


and by solving the equation, yo(t) can be expressed as 


= MeC 
B 4/3Upnor(t + to) i 


Yo(t) (7) 
When we choose 7,3 = 10? as the y at the ICS ini- 
tial point of t = 0, the constant tọ = 1.6 x 107? s. 
Subsequently, the relation between y and ICS compres- 
sion time can be obtained, as shown in Figure 2. From 
the figure, we can observe that regardless of the ini- 
tial y, when the compression time approaches around 
10-4 s, their values are approximately 100. Meanwhile, 
since it is generally believed that y of coherent radia- 
tion in the FRB region is ~ 10t — 10? (Kumar et al. 
2017; Wang et al. 2020), we take its value as Yo,2 ~ 100. 
Furthermore, from Figure 2, we also notice that at a 
compression time of 1074s, the compression length is 
~ 5 x 10?cm. Additionally, it can be found that for 
different initial yrcsi, the compression length remains 
~ 10? cm. Based on these results, we can incorporate 
the values of », » ~ 100, ôt ~ 1074 s, and ôl ~ 5x10? cm 
into the calculations. Because the ôl >> lj, where 
ly ~ 10cm is the bunch size, ôl can be regarded as the 
length of the particle outflow before compression that 
le~ ôl. Here, the physical interpretation of le is that it 
represents the portion of the total particle outflow that 
can be compressed into the size of the bunch. Then we 
can use le to evaluate the volume factor £, 


€ x (Studele)/(Stubel str) ud l./ (cti), (8) 


where lstr is the length of the total particle stream, and 
tiri is the duration of the trigger burst. If the duration of 
the trigger is ti ~ 10 s like neutron star glitch (Ashton 
et al. 2019), then £ ~ 7.7 x 107 19. Through this process, 
a particle stream with ^ 5 x 10? cm can potentially be 
compressed to 10cm in length. Now, we can say that 
the bunch is dynamically formed, and the number of 
particles of each bunch can be estimated as 


Ne ~ 2.5 x 10" B2 15 Vm,13%6 (9) 


Correspondingly, its number density is 
ny = Ne/(YoloStube) ~ 2.5 x 10 Ne 22733 em. (10) 


In the third part, we will calculate the luminosity 
and energy for coherent curvature radiation. For a sin- 
gle electron, the power of curvature radiation (Longair 
2010) is Peer = 2e?yic/(3p?). So, for one bunch, the 
power of coherent curvature radiation is 


2 = 
Pher = Né Peer = NCO sep, (11) 
in which pz ~ 10’ cm is the radius of curvature. If mul- 
tiple bunches (N, ~ 104) are involved in radiation, then 
they are incoherent with each other, which power can 
be written as 


2 
P= 3 Née yoc7 (12) 
Because of the transformation of the coordinates (Ku- 
mar et al. 2017), the isotropic luminosity is Liso = y4 Per 
that 


Liso ~ 3.0 x 107 Np aN? o9% pz? ergs” '. (13) 


Regarding the duration of FRBs, as our model is also 
based on the curvature radiation, some models have al- 
ready analyzed it thoroughly (Wang et al. 2020). There- 
fore, we only provide a brief analysis here that the du- 
ration is related to the height difference (6h) of multiple 
bunches as tpgp ~ óh/c ~ 1ms, and dh is ten times 
the radius of a NS. 


óh = hg — hı € po — pı ~ 10* em, (14) 


where hı (hi) and pı (p2) represent the bunch at posi- 
tion 1 (2), respectively. 

The characteristic frequency (Longair 2010) of the 
curvature radiation is related with y and p that 


v, = 3cy3/ (4np) ~ 0.742 5p7 ! GHz. (15) 


Due to multiple bunches, various y, exist, which will 
move along with magnetic tubes that have different p. 
Hence, the maximum radiation frequency bandwidth for 
one burst is 


3. 3 
v = 3c P2701 Pi Yos. (16) 
Am P1p2 


AVmax xni 


Yo1 and Yo2 are Lorentz factors of bunches at the lowest 
and highest positions within the FRB region (see Figure 
1), respectively. It should be noted that the AYvmas here 
is not the observed bandwidth but rather the theoretical 
bandwidth of a single burst. In reality, only a portion 
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of the bandwidth may be observed due to the signal-to- 
noise ratio. Additionally, the time of bunches emission 
is also sequential, so there may be frequency drift, which 
will discuss later. Above all, the structure of our model 
has been completed, and we will apply it to the FRB 
20190520B in the next Section. 


3. RESULTS AND APPLICATION 


In this section, we will utilize the model derived in 
the preceding section to illustrate its performance under 
various parameters and provide associated inferences. 
Initially, considering our model as a dynamic mecha- 
nism incorporating parameters such as y and p, we will 
present a series of inferences based on the model. Then, 
we will conduct basic simulations to enrich our under- 
standing. Subsequently, these results will be compared 
with the actual observations of FRB 20190520B, with 
the aim of elucidating several phenomena, including but 
not limited to frequency drifting and morphology. 


3.1. Model inferences 


Due to the multiple bunches involved in FRB radi- 
ation, each bunch contributes to the total radiation. 
However, in Eq.13, we assume that their contributions 
are equal (as in most other models). This assumption 
may not reflect physical reality accurately as different 
bunches may originate from varying heights and y. The 
existence of these different parameters is roughly what 
causes the burst duration to broaden to approximately 
~ Ims. As such, it becomes important to examine the 
properties of individual bunches and their contributions 
to the FRB emission. 
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Figure 3. The characteristic frequency distribution diagram 
of a single bunch under different Lorentz factors (yo) and 
curvature radii (p). 


Considering a single bunch, its characteristic fre- 
quency (ve) and luminosity (Lisob) can be computed us- 
ing Ne, Yo, and p according to Eq.11 and Eq.15. It's 
plausible that the higher the launching position of a 


| —— p- yo relation 1048 
1604 
\ 46 
140 10 
120 109 + 
uU 
o D 
> 100 | 1022 S 
$ 
80 E: 
1049 
60 
40 1038 


1.0 1. 1.4 1.6 1.8 2.0 
Curvature Radius (cm) x107 


Figure 4. The luminosity distribution diagram of a burst 
(contains multiple bunches) with different Lorentz factors 
(yo) and curvature radii (p). The blue curve represents the 
relation between p and yo when considering that the value 
of yo will decrease as the p increases. 


bunch (i.e., larger o), the farther it travels, leading to 
longer energy loss and consequently, a smaller yọ value. 
Given these parameters, the variation of ve is graph- 
ically represented in Figure 3. As depicted in Figure 
3, as Yo increases and p decreases, Ve correspondingly 
increases. Concurrently, Figure 4 illustrates the distri- 
bution of Liso across varying values of yọ and p. But, 
as Yo and p increase, Liso also increases, which may ap- 
pear counter-intuitive. This anomaly arises due to the 
fact that both p and yo are provided in the simulation 
code, and the region above the p — Yo curve is physi- 
cally implausible, because when the launching position 
is higher, the yo should be smaller. Hence, only the 
area beneath the curve can accurately reflect the actual 
physical conditions. Having elucidated the properties 
of individual bunches, we can now proceed to scrutinize 
the characteristics of a single burst that encompasses 
multiple bunches. 


3.2. A concise simulation 


In this subsection, we will introduce the concept and 
method of how multiple bunches in different slices are 
simulated. If we assume that the outflow region is cir- 
cular, and the magnetic field lines/tubes are uniformly 
distributed within it, then we can regard that the curva- 
ture radius of adjacent lines/tubes is similar. First, we 
randomly distribute points within a circle of radius Ry 
(outflow region radius) using the Monte Carlo method, 
with the total number of dots equaling the number of 
bunches like N, ~ 10* and each point represent a bunch, 
as shown in Figure 5. Then, we slice the circle uniformly 
into bins along the direction of the line of sight, for ex- 
ample, from right to left, and count the number of points 
(bunches) in each slice. After uniformly dividing the 
emission height, we can calculate several properties of 
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Figure 5. Schematic diagram of the slicing method. Each 
blue point represents a magnetic field line/tube, which can 
also be understood as an outflowing bunch. The area be- 
tween two yellow dividing lines represents a slice, and there 
are a total of 1000 slices in the figure, labeled as S1, S2, 53, 
..., and 81000. The dark green arrow represents the direction 
of the line of sight (LoS), which is the order of the slices. 
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Figure 6. The frequency-luminosity diagram for a single 
burst. Upper panel: dots represent different slices (contains 
multiple bunches), where red-orange points indicate frequen- 
cies higher than 100 MHz and light-orange group meaning 
those below 100 MHz. Bottom panel: The solid blue line 
represents the cumulative luminosity. Note that the vertical 
axis is presented in a logarithmic scale for ease of reading. 


the bunches, such as compression length (Eq.4), coher- 
ent particle number (Eq.9), isotropic luminosity (Eq.13), 
and characteristic frequency (Eq.15). When we assume 
that the Lorentz factor and curvature radiation lumi- 
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nosity of the bunches in one slice is the same in the 
FRB region, we use the sum of the luminosity of all 
bunches to represent the luminosity of a slice (one dot 
in Figure 6). Finally, these luminosities and character- 
istic frequencies are plotted in Figure 6 and divided into 
two groups based on a threshold of 100 MHz. 

For Figure 6, in the upper panel, the luminosity varia- 
tion with frequency is shown for a slice, while the bottom 
panel represents the cumulative luminosity variation for 
one burst. Meanwhile, in our simulation, the high- 
frequency bunches emit radiation earlier compared to 
the middle-frequency and low-frequency bunches. More- 
over, the overall luminosity distribution appears to ex- 
hibit a phenomenon of rapid brightening followed by a 
gradual dimming. If we consider the possibility of a 
low-frequency cutoff and only focus on the middle and 
high-frequency component, the luminosity of a burst 
could potentially follow a power-law distribution. Fur- 
thermore, if the frequency range corresponding to the 
rapid rise in luminosity falls within the observational 
range, our model may potentially explain the occurrence 
of bright spots during the burst. Of course, it should 
be noted that the high, middle, and low frequencies 
are relative values rather than strict bands. Therefore, 
these findings suggest that, for an individual burst, there 
may exist phenomena of frequency downward drifting 
and bright spots. Coincidentally, these features bear 
some resemblance to the observed phenomenon of FRB 
20190520B. 


3.3. Application for FRB 20190520B 


FRB 20190520B, a compelling source discovered by 
FAST (Niu et al. 2022), presents characteristics of ac- 
tive bursts, association with a PRS, and an DMyjyayw that 
markedly deviates from the Macquart relation (Mac- 
quart et al. 2020). These phenomena earmark this 
source as a distinctive and representative FRB, thus 
analyzing this source from a theoretical perspective is 
important. From an energy burst perspective, FRB 
20190520B manifests a log-normal distribution with a 
mean isotropic burst energy of around ~ 10°% erg. As- 
suming a duration of approximately 1ms, the mean 
isotropic luminosity is around 10?! erg s7}, which falls 
within the given value of Eq.13. Considering the higher- 
luminosity bursts within the dataset, which approach 
107? erg s7}, these also conform to the value indicated 
in Eq.13. Concurrently, when considering each slice in 
Figure 6 as having distinct luminosities, the cumulative 
luminosity approximates ~ 1.3 x 109 erg s-!. Limit- 
ing consideration to slices with a frequency above 100 
MHz, the luminosity sum is ~ 8.2 x 104? erg s-1, suf- 
ficient for an FRB 20190520B burst. Even restrict- 
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Figure 7. Dynamic spectra of two bursts from FRB 20190520B. The upper plots mean the pulse intensity varying with time. 
'The bottom-left plots represent the dynamic spectrum after de-dispersion. The bottom-right plots show the intensity varying 
with frequency. Red markers mean frequency channels of the radio frequency interference (RFI), and white horizontal lines 
indicate the areas covered by the RFI. Sub-figure (a) shows the frequency downward drifting of the sub-burst; sub-figure (b) 


shows the bright spot that appears in the mid-frequency. 


ing the calculation to the FAST observation band (1- 
1.5GHz), the cumulative luminosities remain substan- 
tial at ~ 7.0 x 104 erg s-!, thereby still meeting the 
observational requirements of FRB 20190520B. Conse- 
quently, from a luminosity perspective, our model ade- 
quately satisfies observational demands. 

Indeed, only focusing on luminosity or energy analy- 
sis does not entirely illustrate the strength of our model, 
since these aspects can be managed by adjusting model 
parameters. What truly showcases the value of our 
model is its may interpret the complex observational 
phenomena of FRB 20190520B, such as the frequency 
drifting and the localization of bright spots. In regard 
to frequency drifting, we have noted that a downward- 
drifting pattern is frequently observed in FRBs (Pleu- 
nis et al. 2021; Zhou et al. 2022). Wang et al. (2019) 
suggested that for a constant y, this downward drifting 
feature may manifest in FRBs. Our model extends this 
notion by incorporating a variable y, which can also give 
rise to a downward frequency drift phenomenon. In the 
context of FRB 20190520B, we investigated 75 reported 
bursts and identified 12 pairs of sub-bursts with conspic- 
uous downward drifting, akin to the patterns shown in 
sub-figure (a) in Figure 7. If we postulate that different 
sub-bursts from one pair emanate from different compo- 
nents of the same burst, our model could offer a plausible 
explanation for this occurrence. As illustrated in Figure 
6, high-frequency components are emitted earlier than 
their low-frequency counterparts. Further discussion is 
shown in Section 4.3. 


Lastly, we have observed that for some bursts, the lu- 
minosities are not uniform in the time-frequency plot 
but feature bright spots, as exemplified in sub-figure 
(b) in Figure 7. This phenomenon could be ascribed 
to scintillation, which results in conspicuous bright-dim 
stripe patterns in the dynamic spectrum (Main et al. 
2022). However, if bright spots, instead of stripes, are 
observed, they might be associated with the FRB mech- 
anism itself, a common occurrence in FRBs. For in- 
stance, in the case of FRB 20190520B, roughly one-third 
of the bright spots manifest in the middle of the burst 
channel (mid-frequency), as opposed to the top (high- 
frequency) or bottom (low-frequency). In our simula- 
tion, as illustrated in Figure 6, we find that the luminosi- 
ties around 500MHz outshine those in other frequency 
bands. Consequently, this could manifest in observa- 
tional as bright spots occurring near mid-frequencies. 
Meanwhile, for mid-frequencies around 500MHz, the 
CHIME telescope may offer a more suitable platform 
for validating this simulation. Indeed, some bursts 
documented in CHIME/FRB Catalog 1 have exhibited 
a brightening phenomenon at mid-frequencies (Pleunis 
et, al. 2021). It's worth noting that this simulation does 
not precisely match the current observations of FRB 
20190520B. For example, the mid-frequency range with 
the highest luminosity may differ. Nevertheless, when 
examining the region of the FAST frequency in Figure 
6, we can also discern luminosity fluctuations. Thus, the 
luminosity variation in our model may provide a clue for 
explaining the the phenomenon of bright spots. 


4. DISCUSSION 


In this section, we discuss several related implica- 
tions of our model, including model assumptions and 
the initial considerations for bunch formation. Addition- 
ally, possible observational concerns are also addressed, 
which encompass frequency downward drifting, bright 
spots, narrow-band signals, and potential applicability 
to other FRBs. 


4.1. Propagation of FRBs 


A fundamental assumption underlying our model is 
that the radiation emit from the inner magnetosphere. 
This assumption gives rise to a simple but signifi- 
cant question: can the radio signals effectively prop- 
agate through the surrounding medium? If we make 
an assumption that the medium has a density of the 
Goldreich-Julian density (Goldreich & Julian 1969) that 


ngj^10 ?B,Q ~ 10° Bs 153cm’, —— (17) 


where €? is the rotation angular frequency of NS and we 
take Q ~ 1000 rad s7}, then the plasma frequency is 


e?ngj Ms 
Vp,GJ = ( ) ~ 900 GHz. (18) 
TMe 
Moreover, presuming that the density of the surround- 
ing medium is approximately n, ~ 10!" em^?, the cor- 
responding vy, ~ 2.8 x 10? GHz is found. Notably, 
both rj,í; and vp,» are significantly greater than the 
propagation frequency of FRBs (~ 1 GHz). Hence, the 
outflow region must either be located within the region 
of open magnetic field lines, rather than the co-rotating 
region, or a triggering mechanism must be in place to 
clear the surrounding medium and enable the propaga- 
tion of FRB signals. Additionally, several theories pro- 
pose that as long as the condition vprg > V3 [Veye is 
met, radio waves can penetrate high-density medium 
(Arons & Barnard 1986; Kumar et al. 2017), where 
Veye = eBs/(2rmec) ~ 2.8 x 101? B,.45 GHz is the elec- 
tron cyclotron frequency. From this viewpoint, even 
though vrgp < vp,c;, GHz waves can still propa- 
gate that close to the speed of light. Besides, it is 
worth noting that the radiation pressure may exceed 
the plasma pressure in the magnetosphere (Wang et al. 
2022), which could allow FRBs to break out of the sur- 
rounding medium. Thus, based on the mentioned three 
perspectives, the propagation problem of FRB can be 
roughly addressed. 


4.2. Formation of bunches 


Revisiting our initial question that can bunches be 
formed? According to our model, bunches may dy- 
namically form as a result of the electron deceleration 
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mechanism, and the compression effect could result in 
the necessary particle number within the magnetic tube. 
Whether such high-density bunches can form still needs 
further discussion, as many models tend to avoid ad- 
dressing this issue. To begin with, we can provide a 
rough phenomenological analysis from the perspectives 
of electrostatic potential energy and kinetic energy. The 
electrostatic potential energy (Ep) can be viewed as the 
potential energy of all electrons in the bunch with re- 
spect to the electrons entering the bunch, which is writ- 
ten as 

Ep = 3ke? Ne /5r, (19) 


where k is Coulomb constant, and r represents the 
distance between the bunch and the entering electron, 
which is similar to ly. On the other hand, the particle 
kinetic energy (E) refers to the kinetic energy of the 
relativistic electron at the onset of ICS, which can be 
expressed as 


Ey = Eetot — Eo = (mic + peer)? = me, (20) 


where Eetot is total energy of a relativistic electron, 
Eg = mec? is the rest energy of the electron, and 
Pe = ymec is the momentum of the electron. If the 
kinetic energy is greater than the potential energy, it 
can be approximated that the electrons have the ability 
to enter the bunch, thereby facilitating the compression 
process as we described. 

To estimate the balance between the two energies, sta- 
tistical methods can be employed. Assuming that the 
parameters Ne, r, and y conform to a Gaussian dis- 
tribution, a ratio (R) can be introduced to assess this 
matter as 


R= Ex (ty, Oy) / Ep (Uns Hr, On, Or), (21) 


where u and g are mean and variance in the Gaus- 
sian distribution. For this, we can specify values for 
which gy = 1000, 5, = 200, un = 107^, p, = 20cm, 
On = 5 x 1071, and ce, = 10cm. We perform random 
sampling 1000 times for calculation, and the distribution 
is shown in Figure 9. When R > 1, the kinetic energy 
surpasses the potential energy, implying that bunches 
can potentially form. Meanwhile, as Rmean = 1.3 > 1, 
the formation of bunches is feasible in the majority of 
cases. To further mitigate the potential errors caused 
by different samplings each time, we cycle the afore- 
mentioned process 100 times, obtaining the mean value 
of the ratio for each times. Based on these computed 
values, the harmonic mean value is calculated to be 
Rhar = 1.5, closeing with the value depicted in Figure 
9. Hence, from a statistical viewpoint, the formation of 
bunches appears plausible. 
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Figure 8. Dynamic spectra of two bursts from FRB 20190520B. Sub-figure (a) represents narrow-band signal; sub-figure (b) 
means the burst locates at the edge of telescope bandwidth. Other annotations are consistent with those in Figure 7. 
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Figure 9. The histogram of the ratio (R = Ex / Ey) normal- 
ized distribution based on Eq.21. The red and orange dashed 
lines represent the value of R = 1 and R = 1.3, respectively. 


Furthermore, the magnetic pressure generated by the 
strong magnetic fields surrounding the bunches also 
plays a crucial role in confining these bunches. From this 
perspective, we can estimate the relationship between 
the electrostatic potential (Æp) and the work done by 
the magnetic field (Wg). A criterion for when bunches 
may disperse can be determined as follows: 


AE, > AWB, (22) 
where AE, = 3ke?/(5r?) dr and AWg = 4n Par? dr 
represents the changes in electrostatic potential and the 
work done by the magnetic field, respectively. Assuming 
that the bunches are located at an altitude of ~ 107 cm, 
when the surface magnetic field is Bs 15 = 10!? G and 
the neutron star (NS) radius is Ryg ~ 10° cm, the mag- 
netic field within the bunches is Byz = 10!? G, so the 
magnetic pressure Pg ~ B?,/(87). When we substitute 


the parameters mentioned earlier, we find that the cri- 
terion outlined in Eq.22 is not satisfied. Therefore, the 
magnetic field should be sufficient to confine the bunches 
at this altitude. 

According on the above discussions, we can summa- 
rize a few key conclusions. Firstly, an energy-based 
analysis of the formation of bunches under our mech- 
anism provides an intuitive understanding of the pro- 
cess. Secondly, by examining the Eg, the Ep, and the 
We under given parameters, we find that the formation 
of bunches is plausible. Lastly, when we apply statis- 
tical methods to these considerations, the results sup- 
port the likelihood of bunch formation within the de- 
fined parameter space. Additionally, our model does 
not require the bunches to maintain a certain size after 
the FRB emission. Simultaneously, given the exceed- 
ingly short radiation duration of each bunch, there is 
no requirement for them to sustain their structural in- 
tegrity over an extended period (Zhang 2022). However, 
it should be noted that the characteristics of bunches in 
the variable magnetosphere are quite complex, which 
needs further investigation, like considering electromag- 
netic fields, magneto-hydrodynamic waves, and particle- 
in-cell simulation, to make a complete understanding of 
this process. Thus, our discussion represents merely an 
initial exploratory effort concerning this matter. 


4.3. Concerns of the application 


From the standpoint of frequency downward drifting, 
although we have addressed it in Section 3.3, it needs ad- 
ditional consideration. In our model, the prerequisites 
for a downward frequency drift are less stringent, not 
necessitating sub-bursts, and a single burst may also en- 
gender it. This implies that within what is perceived as 


a single burst, multiple structures might exist, although 
they are challenging to distinguish. It should, however, 
be noted that frequency drifting in a single burst might 
also be attributed to fluctuating dispersion measures. 
Especially in the case of FRB 20190520B, with its com- 
plex environmental surroundings, further observational 
and analysis is required. 

'The ability of this model to explain other bursts with 
varying patterns, such as bright spots located in the low 
or high-frequency bands, is something to consider. We 
hypothesize that this could be connected to the struc- 
ture of the outflow region. If we partition the region into 
several equal parts along the line of sight, their horizon- 
tal scales would be the same, but due to the shape of the 
region, their vertical scales would differ. Consequently, 
the area of each slice would also vary, leading to a dif- 
ference in the number of magnetic tubes contained in 
each slice, and this would influence the variation in the 
number of bunches that can be generated by the outflow 
particles. In other words, if the number of bunches in 
a slice is relatively large, the luminosity will be higher 
than in other slices. For a standard circle, assuming a 
line of sight from right to left as in Figure 5, our previous 
analysis suggests that the bright spot should be around 
the mid-frequency band. However, if the region is a tri- 
angle with its apex pointing to the left, the spot should 
appear in the high-frequency band. Conversely, if the 
apex is pointing to the right, the spot should be present 
in the low-frequency band. Hence, based on this con- 
cept, the model could potentially account for different 
burst patterns and even structures with multiple peaks. 

The issue of narrow-band signals is indeed a common 
challenge faced by present relativistic coherent radia- 
tion models. For FRB 20190520B (as seen in sub-figure 
(a) of Figure 8), the majority of observed bursts have 
a frequency bandwidth that is less than half of the in- 
strument bandwidth, i.e., less than approximately 250 
MHz for FAST. In our model, the maximum frequency 
bandwidth is around 3 GHz, which indicates a broad- 
band model. However, a possible explanation is that 
if we assume that the bright part in Figure 6 repre- 
sents the observed bursts with high signal-to-noise ra- 
tios, it could be considered as a narrow-band signal. Fur- 
ther clarification is also required regarding the definition 
of intrinsic narrow-band characteristics, which is neces- 
sary for accurately analyzing the radiation features of 
FRBs. Although certain studies on giant pulses of pul- 
sars and FRBs have suggested that the narrow-band is 
Av/viete ~ 0.1 (Metzger et al. 2019; Thulasiram & Lin 
2021), where Vrele is the bandwidth of the telescope, this 
is still dependent on the instrument and not intrinsic to 
the source. Even different telescope settings and param- 
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eters could affect the identification and classification of 
narrow-band signals (Law et al. 2017). Additionally, we 
find that some bursts of FRB 20190520B are located at 
the upper or lower edges of the telescope’s bandwidth 
as seen in sub-figure (b) of Figure 8. It remains un- 
certain whether there are signals beyond the frequency 
boundary in such cases, and further joint observations 
are needed to provide evidence. If such situations do ex- 
ist, it is important for us to gain a further understanding 
on the morphology and burst mechanism of FRBs. 

Finally, the question arises that can this model be ex- 
tended to other FRBs? We believe that this is currently 
an open question due to the uncertain diversity in the 
origins and mechanisms of FRBs. This model has been 
applied to FRB 20190520B in order to illustrate its ca- 
pabilities in accounting for observed phenomena, but it 
would be a premature claim to consider it as a compre- 
hensive and universally applicable model for all FRBs. 
However, some of the phenomena observed in FRB 
20190520B have also been documented in other FRBs 
(Pleunis et al. 2021), implying possible shared charac- 
teristics among them. For instance, the downward fre- 
quency drifting has been noted in many repeating FRBs, 
including FRB 20121102A (Hessels et al. 2019) and FRB 
20201124A (Zhou et al. 2022). Moreover, for some 
FRBs, the flux of a single burst exhibits non-monotonic 
variation trends with fluctuations or bright spots (Day 
et al. 2020), such as FRB 20180814A (CHIME/FRB 
Collaboration et al. 2019) and FRB 20180916B (Mar- 
cote et al. 2020). Taken as a whole, these observations 
suggest that our model may provide some inspirations 
into explaining these commonly observed phenomena, 
particularly in relation to the compression effects of the 
bunch model. 


5. CONCLUSION 


In this work, we have expanded on the mechanism 
of FRBs, produced by bunches via coherent curvature 
radiation within a NS framework. Our primary focus 
has been on the compression effect exerted by the ICS 
process on outflow particles. We have provided a phe- 
nomenological discussion on the feasibility of bunch for- 
mation from an energy standpoint. Upon further analy- 
sis of the corresponding observational effects, our model 
potentially explains some of the phenomena observed in 
FRB 20190520B. We wrap up our paper with a brief 
summary of key conclusions: 


e We have developed a model for the FRB mecha- 
nism rooted in coherent curvature radiation and 
magnetized NS. This model incorporates the pro- 
cess of bunch compression, thereby offering a novel 
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approach towards understanding the formation of 
these bunches. 


e We introduced the mechanisms of pair production 
and ICS preceding curvature radiation (FRB emis- 
sion). Our estimations suggest that ICS plays a 
significant role in compressing the bunches. Con- 
sequently, these bunches may dynamically form 
and emit curvature radiation along a curved mag- 
netic field tube/family at specific positions. 


e We discussed factors influencing the formation of 
bunches, encompassing the electrostatic potential 
energy within the bunch, the kinetic energy of 
relativistic electrons, and the work done by the 
magnetic field. According to our analysis, the 
compressed bunch model fulfills the conditions re- 
quired for the formation of bunches. 


e We utilized this model to FRB 20190520B and suc- 
cessfully explained various observed phenomena, 


including some basic characteristics and the fre- 
quency downward drifting. Furthermore, we pro- 
vided insights into the position of bright spots in 
FRB 20190520B. 
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APPENDIX 


A. DETERMINATION OF MODEL PARAMETERS 


In this section, we will introduce the selection of pa- 
rameters in the model, as well as the rationale for uti- 
lizing these values. In summary, the parameters can be 
divided into three parts: neutron star (NS) parameters, 
particle outflow parameters, and FRB burst parameters. 
Finally, we will provide a notation list as Table 1. 

Firstly, for the NS part, the parameters include the 
surface magnetic field B,, the volume of magnetic energy 
release Vm, the energy conversion efficiency 7, the ini- 
tial Lorentz factor yg, and the surrounding temperature 
T. Since an FRB-like signal has already been observed 
from Galactic magnetar SGR 1935--2154, it is reason- 
able to consider magnetars as the central powering ob- 
jects. Moreover, the characteristic surface magnetic field 
of magnetars is typically in the range of 1014715 G, and 
considering that FRB bursts are expected to be rela- 
tively violent, selecting a value of 1015 G is justified. 
Regarding the volume part, there are two parameters, 
the radius and depth of energy release that can be trig- 
gered by the mechanism, such as starquakes, magnetic 
reconnection, and explosions. The selection of radius is 
arbitrary, and we assume it may be close to the polar 
cap area of the NS, although they are not necessarily the 
same, which is Re ~ 5 x 10* cm. The selection of depth 
is based on the balance between the magnetic pressure 
and the material pressure that he = B2/(8moncgns) ~ 
4 x 10? em, where one ~ 10!! gem ? is the density of 


NS crust and gns ~ 10!^ ems"? is the gravitational ac- 
celeration of NS. Therefore, the corresponding volume 
is Vm ~ 3 x 10? em?. For the estimation of energy 
efficiency, radio efficiency cannot be directly used. In- 
stead, the efficiency of conversion between matter and 
energy should be selected. Additionally, since the effi- 
ciency corresponding to starquakes and magnetic recon- 
nection is still uncertain, we adopt the efficiency of a 
nuclear explosion 7 ~ 1073. For the initial Lorentz fac- 
tor, we arbitrarily select a value of yẹ ~ 109, which is 
close to the y of particles through the acceleration gap 
(although the meanings of the two are different). Some 
accretion NS systems can produce TeV photons (North 
et al. 1987), which also corresponds to the value of 10°. 
Next is the trigger time, which is selected based on the 
NS glitch time tiri ~ 10s (Ashton et al. 2019). Finally, 
for NS temperature T, it is not the brightness temper- 
ature but the actual temperature of soft photons that 
was heated by the surface of NS and suffuse the inner 
magnetosphere, which we select a value of T ~ 10° K 
(Goldreich & Julian 1969). 

'The second part concerns the outflow parameters, in- 
cluding the ratio of the number of particles from the 
trigger to the outflow region C, the ratio of magnetic 
flux tube cross-section to outflow area &, the Lorentz 
factor after cascading y3, and the compression time ôt. 
The selection of ¢ ~ 107° is flexible because the parti- 
cles generated by the triggering mechanism are isotropic, 
but the outflow region is outside the triggering region. 


So the value is variable, but it should be small. The 
parameter & ~ 10^? is calculated by Stube/Sfiow, where 
Stube = 257 cm? is the cross-sectional area at the bot- 
tom of the magnetic tube and Syiow = 25 x 1087 em? is 
the area of the outflow region. Next is y3 ~ 10?, which 
is related to the pair production and inverse Compton 
scattering (ICS) mechanisms, which will be discussed 
in Appendix B. Finally, regarding the compression time 
ôt ~ 107^ s, we referred to the integrated time of the 
ICS electron lifetime tics ~ YicsMeC?/Peic, where Yics 
we take 73 — Yo,2 ~ 100. 

In the end, we will discuss the parameters related to 
the FRB region, containing the number of bunches Nz 
and the curvature radius p. Here, Ny ~ 3 x 104 is based 
on the average energy of FRB 20190520B, but this pa- 
rameter may vary for other FRBs, so it is also flexible. 
'The selection of the p needs to satisfy two conditions. 
To begin with, it should be larger than the radius of 
the NS (pns ~ 10° cm) because we need the open mag- 
netic tubes in the outflow region. Then, it should be 
larger than the distance that particles move during the 
compression time (~ 3x 10° em). Moreover, considering 
that the compression does not start from the surface of 
the NS, selecting pz ~ 10’ cm is reasonable. 


Table 1. Notation List. 
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Symbol Definition First Appear 
B Magnetic field strength Section 2.2 
Bs Surface magnetic field strength Section 2.2 
c Speed of light Section 2.2 
Eloss Energy loss rate resulting in Section 2.2 
compression 
Eo Rest energy of a electron Section 4.2 
Eetot Total energy of a relativistic Section 4.2 
electron 
Ex Kinetic energy of a relativistic Section 4.2 
electron 
Ep Electrostatic potential energy Section 4.2 
Etot Total energy released Section 2.2 
e Elementary charge Section 2.2 
Ins Gravitational acceleration of NS Appendix A 
he Depth of energy release Appendix A 
h Reduced Planck constant Appendix B 
k Coulomb constant Section 4.2 
Liso Isotropic luminosity Section 2.2 
Lisob Isotropic luminosity of a bunch Section 3.1 
ly Length of a bunch Section 2.2 
le Length that can be compressed Section 2.2 


to la for the particle stream 


'To be Continued 


'Table 1 Continued 
Symbol Definition First Appear 
Is Length of the particle stream Section 2.2 
Me Electron/positron mass Section 2.2 
Neas Cascade number Section 2.2 
Ne Particle number after cascade Section 2.2 
Ntc 'Total outflow particle number Section 2.2 
Niot Initial particle number Section 2.2 
Np Number density in a bunch Section 2.2 
nas Goldreich-Julian (GJ) density ^ Section 4.1 
Pp Pressure of magnetic field Section 4.2 
Pyic ICS power of a bunch Section 2.2 
Pas Power of curvature radiation Section 2.2 
Pric ICS power of a single electron Section 2.2 
Pe Momentum of a electron Section 4.2 
R Ratio between Ej and Ep Section 4.2 
Re Radius of energy release Appendix A 
Rhar Harmonic mean value of Ratio Section 4.2 
Riean Mean value of ratio Section 4.2 
Rys NS radius Section 4.2 
r Distance of electrostatic poten- Section 4.2 
tial energy 
ro Classical radius of the electron Appendix B 
S flow Area of the outflow region Section 2.2 
Stube Cross-sectional area of the mag- Section 2.2 
netic tube 
T NS temperature Section 2.2 
tits ICS electron lifetime Appendix A 
teri Trigger time Appendix A 
Uph Internal energy of the photon Section 2.2 
Vm Volume of the magnetic energy Section 2.2 
release 
v Velocity of the outflow particle Section 2.2 
Qi Interaction angle between the Appendix B 
incident photon and particle 
B Speed ratio between the particle Section 2.2 
and light 
y Lorentz factor Section 2.2 
Yo Lorentz factor of burst Section 2.2 
^ra Lorentz factor after Cascade Section 2.2 
6 Initial Lorentz factor Section 2.2 
AE, Changes of electrostatic Section 4.2 
potential 
AWs Changes of work done by mag- Section 4.2 
netic field 
Av Radiation frequency bandwidth Section 4.3 
Almaz Maximum radiation frequency Section 2.2 
bandwidth for one burst 
oh Height difference of various Section 2.2 
bunches 
él Compression length Section 2.2 


To be Continued 
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Table 1 Continued 
Symbol Definition First Appear 
ot Compression time Section 2.2 
€ Photon energy Appendix B 
¢ Ratio of particle number Section 2.2 
n Ratio of material-energy Section 2.2 
0i Interaction angle between the Appendix B 
scattered photon and incident 
photon 
K Ratio of area Section 2.2 
À Wavelength Section 2.2 
Ing Mean values of Gaussian Section 4.2 
distributions 
Ve Characteristic frequency of cur- Section 2.2 
vature radiation 
Veye Electron cyclotron frequency Section 4.1 
VFRB FRB frequency Section 4.1 
Vi Initial frequency Appendix B 
Vp Plasma frequency Section 4.1 
Vp,b Plasma frequency with n; Section 4.1 
Vp,GJ Plasma frequency with na; Section 4.1 
Vtele Frequency bandwidth of Section 4.3 
telescope 
€ Ratio of volume Section 2.2 
Pus NS radius Appendix A 
One Mass density of NS crust Appendix A 
o Scattering cross-section Appendix B 
Omir Values Variance of Gaussian Section 4.2 
distributions 
Osb Stefan-Boltzmann constant Section 2.2 
OT Thomson scattering cross- Section 2.2 
section 
Q Rotation angular frequency of Section 4.1 
NS 
w Solid angle Appendix B 


Subscript i and o represent the initial and end states of 


the physical parameters, respectively. 


B. CONDITIONS OF ICS 


In our model, we require two effective energy dis- 
sipation processes to diminish the y of the electrons 
prior to the curvature radiation, namely, pair production 
and ICS. Pair production dissipates energy via cascades, 
whereas ICS does so by facilitating interactions between 
high-energy electrons and soft photons. Importantly, as 
y decreases, the mean free path of photons significantly 
increases (Timokhin & Harding 2019), leading to a sub- 


stantial reduction in the efficiency of photon cascades 
during the pair production. Consequently, when ICS is 
dominant, such as when y = 1000, pair production can 
be approximately disregarded. 

Next, we will analyze the conditions of ICS. Consider- 
ing a set of relativistic electrons in a thermal equilibrium 
photon field with a temperature of T' ~ 10° K, the en- 
ergy loss rate due to ICS of the electrons equals this radi- 
ation power, which is Eq. 5. This formula applies when 
the hot photos are not too hard, as yhv < mec?. If not, 
such as y ~ y6, the quantum effects should be consid- 
ered, and the modified equation is (Rybicki & Lightman 
1991) 
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where the €? and € denotes the average of energy squared 
and energy, respectively. At extreme relativity situation, 
the scattering cross-section in unit solid angle should be 
described by Klein-Nishina formula: 


do rae €j € 19 
- o 0 B2 
dw 2e? €o T €i ad , (Bal 


where ro & 2.8 x 1071? em is the classical radius of an 
electron, the €o and e; denotes the energy of final state 
and initial state of photon, respectively, and 0 is the 
angle between the scattered and incident photo. The 
scattering cross-section drops quickly at high frequency, 
and the energy change of each photon is small. Appar- 
ently, the ICS is not an effective energy-losing mecha- 
nism under this condition. However, as a decrease of y, 
the cross-section will increase to ar. When we consider 
that soft photons come from the surface of a NS, then 
the applicable condition of ICS with ør is 


yhw;(1-— coso;) « mec. (B3) 


which gives the restriction relation between the initial 
energy and the initial incident angle (a;) of the photon. 
By considering the photons only coming from the back 
that 7/2 < a; < 7, and the energy of hot photons is 
about 100eV, the condition of y « 5 x 10? is required. 
Meanwhile, in the relativistic electron frame, ICS effi- 
ciency is still low for soft photons coming from directly 
behind the electron, but the angle is relatively small 
~ 1/^, which can be ignored. Thus, it is reasonable for 
us to choose y3 ~ 10? as the end of pair production and 
the beginning of ICS. 
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